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Using a time dependent selfconsistent model for vertical sequential tunneling, we study the appearance of charge instabilities 
that lead to the formation of electric field domains in a weakly coupled doped superlattice in the presence of high magnetic 
fields parallel to the transport direction. The interplay between the high non linearity of the system -coming from the Coulomb 
interaction- and the inter-Landau-level scattering at the domain walls (regions of charge accumulation inside the superlattice) 
gives rise to new unstable negative differential conductance regions and extra stable branches in the sawtooth-like I-V curves. 

PACS numbers: 73.40. Gk, 72.15.Gd 
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Weakly coupled doped semiconductor superlattices are 
an example of non linear systems in which all intrinsic 
properties related to high non linearity such as multi- 
stability, spatio-temporal chaos, etc, can be externally 
modified : e.g. by the application of a bias voltage or 
by the variation of the doping densities in the wells or in 
the contacts. The strong non linear transport that results 
from the Coulomb interaction in these systems presents 
a rich variety of new physical phenomena: multistability 
and electric field domains formation |l|-|3| , self sustained 
oscillations bifurcation to chaos ||, etc. In particu- 
lar, the current flowing through a weakly coupled doped 
semiconductor superlattice presents a complicated saw- 
tooth structure with unstable regions of negative differ- 
ential conductance and multistable regions. This struc- 
ture of branches comes from the charge instabilities that 
appear due to the motion of the domain wall (acummula- 
tion layer) from one well to another. This motion, which 
is due to resonant tunneling between adjacent subbands, 
leads to the formation of electric field domains. In this 
Letter we study theoretically for the first time this phe- 
nomenon in the presence of high magnetic fields parallel 
to the transport direction. New physical questions that 
were not relevant in the absence of magnetic fields can 
be raised. In particular, what happens to the charge in- 
stabilities in the presence of the new energy scale in the 
problem, tuv c ? 

Resonant tunneling through double barriers Jf^dJoC and 
superlattices jyj in the presence of high magnetic fields 
has been widely studied. Also, dynamical instabilities 
and bifurcations to chaos induced by a magnetic field 
in a double barrier resonant tunneling diode have been 
recently predicted fL^ j. It is well known that the ap- 
plication of a magnetic field perpendicular to a two di- 
mensional electron gas produces the formation of Lan- 
dau levels. For ideal samples, the tunneling through 
the heterostructure conserves the Landau level index. 
However, due to interface roughness, impurity scatter- 



ing, or phonon scattering, these conservation rules are 
relaxed and inter Landau level transitions take place 
f|o| P~o|| . Recent magnetotransport experiments on weakly 
coupled doped semiconductor superlattices performed by 
Schmidt et al. ||l3f show new unstable regions of nega- 
tive differential conductance and extra stable branches 
in the I-V curves for a certain range of magnetic fields. 
In this work we will show how the interplay between the 
high non linearity of the system, and inter-Landau-level 
transitions through regions of charge accumulation in the 
structure -domain walls- can explain these new instabil- 
ities. 

The main ingredients of our model are the following: we 
assume that the characteristic time for intersubband re- 
laxation due to scattering which for optical phonon scat- 
tering is about 0.1 ps, is much smaller than the tunnel- 
ing time. The latter is less than 0.5 ns, which is in turn 
much smaller than the dielectric relaxation time respon- 
sible for reaching a steady state, which is about 10 ns for 
the superlattices discussed in Ref. M. This separation 
of time scales, as well as the configuration of a typical 
sample, allows us to assume that only the ground state 
of each well is populated and that the tunneling pro- 
cesses are stationary. These assumptions then justify the 
use of rate equations for the electron densities at each 
well. In order to include relaxation in the transport di- 
rection due to the different scattering mechanisms, we 
suppose that the spectral densities in the quantum wells 
are Lorentzians whose half width 7 is a parameter re- 
lated to the scattering life time, which is of the order of 
picoseconds. Relaxation in the planes perpendicular to 
the transport is considered imposing current conserva- 
tion in the stationary limit. This condition determines 
the sequential current as well as the Fermi energies within 
each well {e Wi } |lj,|l^]. The electrons are described by a 
Fermi-Dirac distribution in each well, which implicitly as- 
sumes that some scattering mechanism thermalizes them 
towards a local equilibrium situation. In the absence of 



such a mechanism for thermalizing the electrons the dis- 
tribution functions are nonequilibrium objects and the 
rate equation is equivalent to a quantum kinetic method 
-i.e, Keldysh or Kadanoff-Baym, for example. The rate 
equations -for a given set {e Wi }- are: 



dm 

at 



i = 1, . . . ,N. (1) 



In this equation, N is the number of wells, Jj,t+i are 
the interwell currents and Jo,i and Jjv,jv+i are the cur- 
rents through the contacts (emitter and collector re- 
gions). These currents are calculated by means of the 
Tunneling Hamiltonian method p3, and have depen- 



dencies J, 



M+i 



Ji. 



M+ i(e Wi ,e Wi+1 ,{$}). {$} denotes the 
set of variables coming from the electrostatics: potential 
drops at the accumulation and depletion layers, barriers 
and wells, etc. These must be calculated self consistently 
for each applied bias. More details of the selfconsistent 
procedure and of the electrostatic model considered here 
are given in Ref. [ fL5[ . These rate equations, which in- 
clude all of the relevant dynamics of the selfconsistent 
problem, can be rewriten in terms of an Ampere's law 
that explicitly shows the time-dependent current con- 
sisting of a displacement current plus a tunneling term 

© lit + J *-M = J W- Here T is the electric neld 
at the i-th barrier and e is the GaAs static permittiv- 
ity. Nonetheless, in this work we restrict ourselves to the 
stationary regime — ► 0. In this limit, the interwell 
currents and the currents from the contacts are all equal 
to the total sequential current J. Following Ref. and 
considering the formation of Landau levels in the planes 
perpendicular to the transport, we obtain the following 
expressions for the tunneling currents: 
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In these expressions f £u . (e z ) = (1 + e k B T ) 1 are 
the Fermi functions, e(n) = huj c (n + h), w c = ^| is 
the cyclotron frequency, n runs over all Landau levels 
below the Fermi energy of the region (emitter or i-th 
well) from which the electrons are tunneling. The Fermi 



energies, e Wj , must be calculated self consistently. Also, 
g = is the factor of degeneracy for a given magnetic 
field B and area S, j labels the resonant state in each 
well i with energy e % C p m is the total number of reso- 
nant states within the well contributing to the current 
and Bi,i + i = ki/(vu + a^ 1 + Q^i), where ki and are 
the wave vectors in the wells and the barriers, respec- 
tively. These depend on the local electric field. Finally, 
Ti = 16/c i _ifc l afe~ 2Ql!i (A;2_ 1 + a?)" 1 (A;? + a})- 1 is the 
dimensionless transmission probability through the ith 
barrier, and w and d are the well and barrier widths re- 
spectively. 




FIG. 1. I-V curves for different magnetic fields in in- 
crements of 0.5T. The curves are shifted in intervals of 



0.2 x 1(T 



for clarity. 



As we mentioned above, we use Lorentzians for the 
spectral densities in the wells A l C j(e) — 7/[(e— £cy) 2 +7 2 ]- 
Although this assumption can be improved by perform- 
ing a calculation of the self-energies including microscop- 
ically the scattering (i.e., including the energy and mag- 
netic field dependence of the selfenergy), previous micro- 
scopic calculations of impurity scattering in the presence 
of magnetic fields have shown this Lorentzian form to be 
a good approximation J^] . On the other hand, our hamil- 
tonian restricts inter Landau level scattering just within 
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the wells. The effect of this contribution on the current 
has been shown to be the most important one of the pos- 
sible impurity scattering || and interface roughness ||] 
contributions in a double barrier. 

Thus, our phenomenological model for the scattering ac- 
counts for the main physics: inter-Landau-level scatter- 
ing at the domain wall (see below), and makes the non 
linear problem tractable numerically. Also it is impor- 
tant to emphasize that the Fermi energies in each well 
are variables; i.e, after each tunneling event the lateral 
energy e(n) is not necessarily conserved. After reach- 
ing a stationary state, at a certain bias, some amount 
of charge is accumulated in the well where the domain 
wall is formed. Depending on the degeneracy, it could 
imply an increase in the maximum Landau level occu- 
pied in one well with respect to its neighboring. The 
interwell tunneling then could involve processes where 

An n max ~ n max ¥ z ^- 

We have analyzed a superlattice consisting of 15 wells 
of GaAs of 90 A thickness and GaAlAs barriers 50 A 



ide. The emitter doping is No = 2 x 10 



and 



the doping in the wells is = 1.5 x 10 11 cm -2 , with 
7 = AmeV , and T = OK. For zero temperature the ex- 
pression for the charge density in the i-th well is analytic: 

-e(n) 
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In Fig. 1 the I-V curve is plotted for different magnetic 
fields (the curves are shifted in the abcissa axis for clar- 
ity). For low magnetic fields, we observe a main peak 
at low bias coming from ground to ground state inter- 
well tunneling, and a sawtooth structure at high bias 
due to electric field domains formation. As the magnetic 
field increases, additional unstable negative differential 
conductance and stable branches show up in the current 
density. 




FIG. 2. Differential conductivity versus bias voltage and 
magnetic field. 

This extra structure dissappears as the magnetic field 
is increased above ~ 13T. In order to study the range 
where this effect manifests itself in more detail , we plot 



in Fig. 2 a contour plot of the differential conductiv- 
ity as a function of B and the external bias V (negative 
differential conductance regions are the white lines). 
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FIG. 3. a) Electric fieldTTor Vi = 0.451Volt and 
Vt = 0.496Volt, b) Electric fields for V a = 0.470Volt and 
V b = 0.475Volt, c) Detail of the I-V curve for B=11T show- 
ing the extra branch. 

One can see that the additional structure in the current 
takes place between 8 and 12.5 T where the differential 
conductivity shows two branches at fixed magnetic field 
that repeat periodically. This result is in good qualita- 
tive agreement with experiment, that shows this effect 
up to 19 T [Q. At lower magnetic fields around 6.5 
Tesla, and voltages close to 0.4 Volt, there is a small 
structure reflecting the formation of a multistable solu- 
tion (see Fig. 2). Our calculations have been performed 
by direct numerical integration of the dynamical system 
|15| in order to permit a direct comparison with the ex- 
periment [jl3|]. The experimental I-V curves have been 
obtained by sweeping the DC bias only in one direction 
-i.e, by increasing the DC bias in between the emitter 
and the collector. Nonetheless, it is possible to use more 
complicated numerical techniques such as the numerical 
continuation method in order to obtain all of the un- 
stable and multistable regions. The multistable regions 
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for the current could be observed in the experimental 
curves along up and down voltage sweeps. 
We shall concentrate in one region of bias for the current 
density and fix B=11T in order to obtain an explana- 
tion for these new structures appearing at intermediate 
magnetic fields. We plot in Fig. 3c, a blow up of the 
current for B=11T in the region of voltages where one 
third of the wells already belongs to the high field do- 
main. We have analyzed the electric field distribution 
for two different biases that give the same current den- 
sity: Vi = 0.451Volt and V 2 = 0.496Volt (see Fig. 3a). 
The first voltage Vi, corresponds to a current solution 
that belongs to the first branch and the second one V2, 
corresponds to a solution in the second branch (Fig. 3c). 
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FIG. 4. Detail of the calculated electrostatic profile for 
0.495 Volt and B=11T. Dotted lines are the Landau levels. 
Also the calculated maximum occupied Landau levels for each 
well are shown. 

As we can see, the current density presents a disconti- 
nuity at V a = 0.47Volt, and it drops abruptly. Increasing 
the voltage, an additional branch shows up. The current 
decreases from one branch to the next due to the electric 
field distribution (see in Fig. 3b the electric field dis- 
tribution in the neighborhood of the current discontinu- 
ity). For voltages higher than V a = 0.47Volt the resonant 
states in neighbor wells at the domain wall become dis- 
aligned in energy, their wavefunctions overlap decreases, 



and the current decreases as well. The current disconti- 
nuity between both branches is due to the strongly non 
linear charge and field distributions. The electric field is 
discontinuous between V a = 0.47Volt and Vb = 0.475Volt 
at the 9th barrier (Fig. 3b). A detail of the calculated 
potential profile for a voltage slightly less than V2 (Fig. 4) 
shows the electrostatic profile at the domain wall region. 
This voltage corresponds to the formation of the extra 
branch. The position of the resonant states -which have 
been drawn as discrete ones for simplicity- and Landau 
levels energies at the domain wall are depicted as con- 
tinuos and dotted lines respectively. Also, the calculated 
maximum Landau level partially occupied in each well, 
n max = anQl n Li=l respectively are indicated. These 
are obtained from the calculated Fermi levels for this volt- 
age. We conclude that the extra branch of the current 
is a result of inter Landau level transition at the domain 
wall , An = v? max — v? max = 1, between the 8th and 9th 
wells. In this situation, the electrons tunnel from a re- 
gion with low charge density accumulation towards the 
domain wall while changing their maximum parallel en- 
ergy to Twj c (n 9 max + \). As the DC voltage increases, the 
domain wall moves from one well to the previous one -in 
this particular case from the 9th well to the 8th well- 
and the structure discussed above repeats periodically 
(see Fig.I). For magnetic fields below 8 Tesla, the mo- 
tion of the domain wall does not involve a change in the 
maximum occupied Landau level and there is no inter 
Landau level scattering involved. In the regime of high 
magnetic fields, where the domain boundary enters the 
magnetic quantum limit, the lowest Landau level is the 
only one occupied in the wells because high degeneracy, 
and we obtain only one branch which repeats periodically 
as the domain wall moves. 

The reason why we do not obtain exactly the same range 
of magnetic fields as in experiment [[[3) , is mainly due to 
the simplicity of our model for scattering as we have dis- 
cussed above. We have observed that by increasing the 
doping in the wells, the range of magnetic fields becomes 
closer to the experimental one. Furthermore, we have not 
included spin-splitting and exchange effects in the calcu- 
lation. At high fields, the exchange could be important 
and the spin splitting must be considered. This effect 
would also give new features in the non linear current at 
high fields. This will be addressed in future work. 
In summary, we have proposed and solved a dynamical 
selfconsistent model for studying the sequential transport 
through a doped superlattice in the presence of magnetic 
fields. We study the appearance of charge instabilities 
leading to the formation of electric field domains. For in- 
termediate magnetic fields, the appearance of new unsta- 
ble negative differential conductance regions and new sta- 
ble branches in the sawtooth-like I-V curves is discussed. 
These new features in the current can be explained in 
terms of inter Landau level scattering occurring at the 
domain walls. 
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